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Abstract 

The problem of creating well-collimated beams of atoms escaping from 
a trap is studied. This problem is of high importance for the realization of 
atom lasers. Nonadiabatic dynamics of neutral atoms in nonuniform magnetic 
fields, typical of quadrupole magnetic traps, has been considered. The main 
result of this report is that an unusual semiconfining regime of motion is 
found, when atoms are confined from one side of an axis but are not confined 
from another side. This regime can be achieved by means of only magnetic 
fields. The formed atomic beam can be forwarded in arbitrary direction. 
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1 Introduction 



The possibility of realizing Bose-Einstein condensation in trapped dilute gases (see 
reviews [1,2]) demonstrated that a macroscopic number of bosons could be produced 
in a single quantum state of trapped atoms. The occupation of a single quantum 
state by a large number of bosons is the matter-wave analog of the storage of photons 
in a single mode of a laser cavity. The device that could emit coherent beams of Bose 
atoms, similarly to the emission of photon rays by light lasers, has been called atom 
laser [3-10]. Output couplers that are used for extracting condensed atoms from a 
trap are based on transferring atoms from a magnetically trapped state to another 
internal state that is not trapped [11-14]. The transferring field can be either a 
weak radiofrequency field including spin flips between trapped and untrapped states 
[11-14] or a field of a laser beam stimulating Raman transitions between magnetic 
sublevels [15]. In the first case, atoms in the output state simply fall down from the 
trap under the action of gravity. In the second case, when the atoms are transferred 
to a magnetic untrapped state, they are also given a momentum kick from the 
photon recoil, so that they exit the trap in a well-defined beam whose direction and 
velocity are governed by the irradiating Raman laser. 

Since the very first condition on any laser is that its output be highly direc- 
tional [9], it is important to develop mechanisms permitting to create well-collimated 
beams of atoms coupled out of a trap. A novel mechanism of creating highly direc- 
tional beams if atoms leaving a trap has recently been suggested [16] and different 
regimes of operations have been studied [17-20]. 

This mechanism is based on transferring trapped atoms to a nonadiabatic initial 
state, which can be done by a short triggering pulse. There exist such initial spin 
states for which the motion of atoms becomes semiconfined for the same configura- 
tion of the trap magnetic field. Then atoms quit the trap flying out predominantly 
in one direction and forming a well-collimated beam. In this report, we give a brief 
outline of the semiconfining regime of motion. 



2 Solution of Differential Equations 

The main equations of motion for an atom of mass m and magnetic moment fiQ 
are usually written in the semiclassical approximation for the quantum-mechanical 
average of the real-space coordinate, R= {R a }, and for the average S— {S a } of the 
spin operator, with a — x,y, z. The first equation writes 

d 2 R a _ /fo d B . s 

dt 2 ~ m b 'dR a ' 1 j 

where B is a magnetic field and, for a while, we omit the gravitational force and a 
collision term whose role will be discussed later on. The equation for the average 



2 



spin is 

d S _ Ho 

dt h 



S x B ■ (2) 



The total magnetic field B—Bi + B 2 consists of two terms. The first is the 
quadrupole field 

Bi= B[ (R x e x +R y e y +XR Z e 2 ) , (3) 

which we write in a slightly more general form than in Ref. [16], including the 
anisotropy parameter A. The latter was taken in Ref. [16] in the standard form as 
A = —2. But, if several magnetic coils are involved in the formation of the field (3), 
the anisotropy parameter A can be varied. The possibility of such a variation will be 
important in what follows. The second part of the magnetic field B is a transverse 
field 

B 2 = B 2 h (t) , h (t) = h x e x +h y e y , (4) 
in which h a = h a (t) and | h (t)\ — 1. In Ref. [16] the rotating transverse field, with 

h x = cos out , h y = sin cut , (5) 

was considered, although, as will be shown here, this is not the sole possibility. 
Introducing the dimensionless space variable r=R / Rq = {x, y,z}, measured in 
units of the length Rq = B 2 /B[, and defining the characteristic frequencies 



" 1 = V^ " 2 = ^T' (6) 



uil (S x ) (7) 



we may rewrite Eq. (1) as 

d 2 7 

dt 2 

and the spin equation (2) as 

^ = MS, (8) 

where the matrix A = [A a p] is antisymmetric, A a p = —Ap a , A aa = 0, and 

A 12 = Xz , A 23 = x + h x , A 31 = y + h y . (9) 

Keeping in mind the following discussion on the generality of the approach, we do 
not require here that h x and h y in Eq. (9) be necessarily of the rotating form (5). 
What is required is the existence of small parameters, one of which is 

^ « 1 (10) 

UJ 2 
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and another one is 

u ^ i - d h 

— <C 1 , to = max 

u 2 1 



dt 



(11) 



In the particular case of the rotating field (5), the value of u in Eq. (11) coincides 
with the rotation frequency. 

It is worth noting that the inequalities (10) and (11) are assumed from the very 
beginning. Numerical estimates of Ref. [16] give uoi/uo 2 ~ 1CT 5 and uo/uo 2 ~ 1CT 3 . 
The existence of the small parameters (10) and (11) allows us to apply for solving 
the system of equations (7) and (8) averaging methods [21-23] as developed for 
multifrequency systems [24,25]. A generalization of averaging methods, called the 
scale separation approach [26], has been developed for treating statistical systems 
that can be described by stochastic differential equations. However, while equations 

(7) and (8) contain no stochastic terms, the mathematical foundation of solving 
them is directly based on the averaging methods [21-25]. 

According to the small parameters (10) and (11), the functional variables r 
and h are to be treated as slow compared to the fast variable S- This means 
that r and h change very little, being almost constant, during the oscillation time 
of S- Therefore, r and h can be regarded as quasi-invariants with respect to 
S- Such quasi-invariants in the theory of classical Hamiltonian systems are called 
adiabatic invariants [25] , whose existence should not be confused with the adiabatic 
approximation. Following the averaging methods [21-25], we can solve equation 

(8) for the fast spin variable with the slow variables r and h treated as quasi- 
invariants. The difference between the latter is that r , as a function of time, is 
defined implicitly by Eq. (7), while h is an explicitly given function. Of course, this 
difference is merely technical but not principal, since both r and h, according to 
inequalities (10) and (11), are slow variables. 

The matrix A defined in Eq. (9) depends on time only through r and h- This 
permits us, introducing the small parameter 



(12) 

I UJ 2 UJ 2 ) 



e = sup 

to write the evolution equation 

in which C is a matrix with bounded elements for any e < 1 and all t > 0. Conse- 
quently, the matrix A has also to be treated as a slow matrix variable. This allows 
us to look for a particular solution of the spin equation (8) in the form 

Si (t) =b l (t)exp{^(t)} , (14) 

in which bi, with i = 1,2,3, are the eigenvectors of the matrix A, given by the 
eigenproblem A bi— oti bi- The nice antisymmetric structure of A makes it possible 
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to easily organize that the set {bi} would form an orthonormal basis. Substituting 
expression (14) into Eq. (8), we find 

(fi(t) = (u 2 ai- bi - bij dt , (15) 

as in Ref. [1]. As is clear, the form (14) is not, strictly speaking, an exact solution 
of Eq. (8), but this is what in the theory of differential equations called an asymp- 
totically exact solution with respect to the small parameter e, in the sense that this 
solution becomes exact when e — > 0. The latter is evident from Eq. (13), since, 
when e — > 0, the matrix A becomes a matrix with constant coefficients. 

In the guiding center approach [27], which is a variant of the averaging methods 
[21-25], the solution of type (14) is called the guiding center because it describes 
the leading behaviour of an exact solution. Corrections to the guiding center can 
be obtained by presenting the exact solution as an expansion around this guiding 
center. Such expansions, because the guiding center itself is generally a function of 
the small parameter, are called the generalized asymptotic expansions. These have 
been introduced and studied by Lindstedt and Poincare [28] and are actively used in 
averaging methods [24-26]. Corrections to the guiding center describe small fastly 
oscillating ripples of higher harmonics and with the amplitudes of the orders of in- 
creasing power of the small parameter. As is well known [29-31], the construction 
of generalized asymptotic series is not unique, being based on an incomplete ex- 
pansion in powers of the small parameter. For instance, one can separately expand 
amplitudes and phases, as in the Lindstedt-Poincare technique [28-31]. In such in- 
complete expansions, one may keep different number of terms in the amplitudes and 
phases, which is dictated by the properties of the particular equation. For example, 
this can be necessitated by cancelling secular terms or by considering more accu- 
rately the phase that, entering in an exponential, can influence the solution stronger 
than corrections to the amplitude. In our case, the nonuniqueness of defining the 
guiding center (14) could result in the possibility of writing, instead of the phase 
(15), the phase 

(fi(t) = uj 2 a.i{t)dt . (16) 
J o 

This is admissible, since the second term in the integrand of Eq. (15) is small as 
compared to the first one. The sole doubt could be whether it is possible to omit 
the second term in Eq. (15) if one of the eigenvalues ctj is zero. In our case it is 

really so for a 3 = 0. However, again owing to the nice antisymmetry property of A, 

— »* — » — >* — > 

we can show that b\—b2 and 63=63. Then, because of the normalization condition 

I ^3 I 2 = bz= 1, we have 2 6 3 -& 3 = d\ 63 \ 2 /dt = 0. Thus, y2 3 (t) = exactly. Which 
of the expressions, (15) or (16), is to be used makes no difference for what follows, 
since at the next step of the method [24-26] one has to substitute the solution for 
the fast variable into the right-hand side of the equation for the slow variable and 
to average over time this right-hand side. In the course of this averaging, all fastly 
oscillating terms disappear, irrespectively to the fine structure of their expressions. 
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When speaking about asymptotically exact solutions, one should not confuse 
this with the asymptotic behaviour with respect to small time t — > 0. The solutions 
we are talking about are asymptotic with respect to small parameters, such as (10), 
(11), or (12). These parameters are always small for any t > 0. There is also the 
following invariance property of Eq. (8) if A is an antisymmetric matrix. Let [SV/(t)] 
be the fundamental matrix of solutions to Eq. (8). Then [32], one has 

det[Sij(f)] = det[Sij(* )] exp |cu 2 J* Tr A(t')dt' 

for any to- For an antisymmetric matrix, TrA = 0. Hence, det[SV,(t)] = const. 

Thus, if [S'jj(to)] is the fundamental matrix at some t — to, then the solutions Si (t) 
are linearly independent for any t > 0. The total solution of Eq. (8) is the linear 
combination 

3 

s (t) = J2 a iSi (0 ( 17 ) 
i=i 

of linear independent particular solutions, with the coefficients a, —So • hi (0) de- 
fined by initial conditions. In our case, the asymptotically exact particular solutions 
are given by Eq. (14). 

Let us emphasize that all the consideration is based on and justified by the aver- 
aging methods [21-26] requiring the existence of small parameters that permit one 
to classify functional variables onto fast and slow. Without such small parameters 
we would have quite a different story. 

After the solution (17) for the fast variable is found, then following the multi- 
frequency averaging methods [24-26], one has to substitute it into the equation (7) 
for the slow variable and to average over time this right-hand side, which gives 

d 2 7 - -> - - 

=F= oof < S x e x +S y e y +\S Z e z > . (18) 

Till now, the concrete form of the transverse field (4) has been of no importance, 
provided that condition (11) holds true. To explicitly accomplish the time averaging 

in Eq. (18), we need to specify the field h- With the rotating field (5), we find 
F= \hu\ [(l + x)S° x + yS° y + \zS°] (x e x +y 7 y +2\ 2 z e 2 ) , 



where 



fi = (l + 2x + x 2 + y 2 + AV)(1 + x 2 + y 2 + \ 2 z 2 ) " J ' 



As is shown in paper [16], the semiconfining regime of motion can be realized 
taking the initial spin polarization given by the initial conditions 

Si = S° y = , S° z = S ? . (19) 
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Such initial conditions can be prepared in several ways. The first possibility is to 
confine atoms in a trap, in which all atoms are polarized having their spins in the 
z direction, which e.g. can be done by means of the trap of Ref. [33], being a 
quadrupole trap with a bias field along the z axis. Then the longitudinal bias field 
is quickly switched off and, at the same time, a transverse field, as in Ref. [34], 
is switched on. Thus, we get the desired initial conditions. Another possibility is 
to prepare spin polarized atoms in one trap and quickly load them into another 
trap with the required field configuration. The possibility of realizing two ways of 
transferring atoms from one trap to another, by means of sudden transfer as opposed 
to adiabatic transfer, is discussed in Ref. [35]. The third possibility could be by 
acting on the trapped atoms with a short pulse of strong magnetic field, polarizing 
atomic spins in the necessary way. 

With the initial conditions (19), equation (18) reduces to the system 

^ = \su\hzx , ^ = \*Sulf lZ * , (20) 

the equation for the y component having the same form as for the x component. The 
system of equations (20) can be solved analytically only for the case | r | 1 when 
fi — > 1, as is considered in Ref. [16]. This consideration proves the appearance of 
the semiconfining regime of motion and the existence of collimation with the aspect 
ratio 

^~l^o| 3/2 , 
z(t) 

as t — > to, where to is defined by the initial conditions for the space variable. When 
| r | becomes not small, one has to resort to numerical solution which has been 
done in Refs. [18,20] confirming the existence of the semiconfining regime and the 
formation of well-collimated atomic beams. 

Let us stress again that asymptotic solutions obtained by the averaging methods 
are asymptotic with respect to small parameters and have nothing to do with short- 
time approximation. As follows from the theory of the averaging methods [21-25], 
the latter provide accurate solutions to differential equations for long time intervals. 



3 Applicability to Different Fields 

The semiconfining regime of motion can be realized for a wide class of magnetic 
fields. Not only the anisotropy parameter A in the quadrupole field (3) can be 
varied, but also different transverse fields (4) can be employed. In the previous 
section, the particular form of a transverse field has not been specified till formula 
(18), showing by this that the same consideration is valid for any transverse field 
satisfying condition (11). Now we shall show that the choice of the rotating field (5) 
is not compulsory and other transverse fields can be taken provided condition (11) 
holds true. 
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Let us consider the simple case of a constant transverse field 



h x = const , h y = const , (21) 

where, similarly to Eq. (4), it is assumed that | h | = 1- Then the small parameter 
(11) is identically zero. 

With the transverse field (21), the right-hand side of Eq. (18) becomes 

F= hu\ \{x + h x )Sl + (y + hy)S y ) + XzS°] [(x + h x ) e x +(y + h y ) e y +X 2 z e. 

(22) 

where 

f 2 = [( x + h x ) 2 + {y + h y ) 2 + X 2 z 2 ]~ 1 . 
Accepting the initial conditions (19), we come to the equations 

^ = XSu 2 J 2 z(x + h x ) , ^ = X 3 Su 2 J 2 z 2 . (23) 

One may notice that these equations are invariant with respect to the change A — > 
—A, S — > —S or to the change S — > —S, z — > —z. Therefore, it would be sufficient 
to consider a case of one fixed sign for XS, say, XS > 0. Passing to another sign of 
AS* < is equivalent to the inversion z — > —z. 

Equations (23) are very similar to Eqs. (20). In the same way as in Ref. [16], 
we can demonstrate that the solution to Eqs. (23) corresponds to the semiconfining 
regime of motion. For | r \ <^ 1, when f% — > 1, we have 

^ = XSu\h x z , ^ = X 3 Suj 2 z 2 . (24) 

The second equation here is the same as in Eq. (20) if f± — > 1. The analysis shows 
that the motion along the z axis is semiconfined. The collimation of the forming 
directed beam is characterized by the aspect ratio 



z(t) 



~ \t-t \ 2 \n\t-t \ , (25) 



as t — > to- Comparing this with the aspect ratio for the rotating field, we take into 
account that for r = \t — t \ 1 one has t 1 / 2 \ lnr| 1. Therefore, the constant 
transverse field provides even better collimation than the rotating field. 



4 Influence of Gravitational Force 

Atoms having mass are certainly subject to the action of the gravitational force. For 
the adiabatic motion of atoms confined in a trap, this force is not as significant since 
it can always be easily compensated by an additional gradient magnetic field. But 
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what is the role of this force in a nonadiabatic motion? The principal answer to this 
question is, actually, almost evident: If the gravitational force can be compensated 
by additional gradient fields, it will lead to no principal changes in the motion of 
atoms inside a trap. However, to be precise, let us turn to mathematics. 

Recall first of all that the coordinate axes everywhere in our formulas have been 
related to the magnetic field configuration, so that the axis z is an axis of the device, 
but not necessary the vertical axis related to gravity. As far as the device can be 
oriented arbitrarily, the gravitational force can also be directed along different axes. 
For instance, assume that this force is — mg e x , where g ~ 10 3 cm s~ 2 is the 
standard gravitational acceleration. Adding this force to the right-hand side of Eq. 
(18) results, instead of Eqs. (20), in the equations 

^ = a;? [±Shzx - a) , ^ = \*Sulh* , (26) 
in which the notation 

a =iti (27) 

is used. The occurrence of an additional term in the first of equations (26) will lead 
to the gravitational drift along the x axis. To compensate the gravitational force, 
an auxiliary gradient magnetic field is to be superimposed. This is equivalent to 
the increase of the coefficient near e x in the quadrupole field (3), which yields the 
scaling x — > \ x x in Eq. (26). One may also increase the anisotropy coefficient A. The 
final aim is to make the ratio a/\\ x small, which thus reduces the influence of the 
gravity force. This ratio is actually small already for existing traps, without involving 
additional compensating fields. For instance, taking the values iO\ ~ (10 2 — 10 3 )s _1 , 
Ro ~ (0.1 — 0.5)cm, |A| = 2, and A^ = 1, typical of many modern traps, we have 
a/\\ x of the order 10 _1 — 10~ 3 . Increasing the parameters A and \ x , by switching 
on additional gradient fields, the ratio a/\\ x can be made arbitrarily small, thus 
making the influence of the gravitational force negligible. 

Moreover, we have an additional possibility of changing the orientation of the 
device with respect to the gravitational force. For instance, the latter can be directed 
along the z axis, so that this force is — mg e z . Then, instead of Eqs. (26), we get 

^ = ^Sulhxz , ^| = ul(\ 3 Sf lZ 2 - a) . (28) 
dt 2 2 lJ ' dt 2 1V J ' K J 

From here it is seen that the influence of the gravitational force can be reduced by 
increasing the anisotropy parameter A. To demonstrate this more accurately, let us 
consider the same case as earlier, when | r \ <^ 1. Then, integrating once the second 
equation in the system (28), we have 



dt 3 1 



z 3 -4-^j(z-Zo)+( 
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where 



3i 2 



The latter equation, by means of the substitution 
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z(t) = ^V(r -r ), t = U!t , (29) 



can be transformed to the Weierstrass equation 




2 

= 4P 3 - g 2 V - g 3 , 



with the Weierstrass invariants 

.. _ (A 3 ^) 3 / 3 3n 



#2 



From the properties of the Weierstrass function V(t) it follows that the semiconfining 
regime of motion is realized when g\ < 27g%, which gives 

4« 3 / q 3a A 2 

< N-T^o-C • (30) 



A 9 S 3 V u X 3 S 

As is evident form this inequality, the terms related to the gravity can be strongly 
reduced by increasing the anisotropy parameter A, that is by increasing the field 
gradient in the z direction. Even for the typical value |A| = 2, without involving 
additional compensating fields, the left-hand side of the inequality (30) is of the 
order 10~ 3 a 3 . For a ~ 10~ 3 this yields 10~ 12 , which is very small. Increasing the 
parameter A can practically completely eliminate the influence of gravity. Note also 
that for the case AS* < 0, the inequality (30) is always valid since its left-hand side 
is negative while the right-hand side is nonnegative. The latter case shows that the 
gravitational force even may help to realize semi confinement, which happens when 
the direction of this force coincides with that of the atomic escape. In this way, it 
is not difficult to prepare such magnetic fields and to choose the orientation of the 
device so that the gravitational force would not essentially disturb the semiconfining 
regime of motion. 



5 Role of Random Collisions 

In order to investigate the role of atomic collisions, we have to compliment the 
right-hand side of the evolution equation (18) by an additional term describing 
these collisions, so that we can write 

d 2 r^ -» 

-ST =F +7 { , (31) 
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where F is the same expression as in Eq. (18), 7 is a collision rate, and £ is 
a vector whose properties are defined by the details of the pair collisions. Two 
principally different situations can exist. One is when the effective collision force 

7 £ is larger than or comparable with F- Then it is clear that the atomic motion 
is essentially governed by the collision force whose particular properties become 
important. However, this case is of no interest for us since, as is evident, in the 
presence of strong and frequent collisions no directed semiconfined motion can occur. 
Any organized directed motion will be strongly suppressed by disorganizing random 

collisions. The second situation is when 7 £ is much weaker than F- And solely this 
case is of interest since only then an organized unidirectional motion can arise. But 

when the force 7 £ representing random pair collisions is just a weak perturbation, 
then the details of this force are not of great importance and it can be modelled by 
the stochastic white-noise variable. It is this approach that was employed in Ref. 

[16], where some additional simplifications for £ were assumed. These simplifications 
are not principal, and the consideration can be easily generalized to the case of an 

anisotropic random vector £ = where fx = x,y,z. In the general anisotropic 

case, the stochastic properties of the set {£^} of white-noise random variables is 
characterized by the stochastic averages 

»=0 (fi = x,y,z), <^(*)^(0»=2^M(*-0, ( 32 ) 

in which is a diffusion rate in the fx direction. 

Adding the anisotropic random force to Eq. (20), for | r | <g; 1, we have 

^ = ^Su 2 1 zx + 1 £ e , ^ = X 3 Suj 2 1 z 2 + 1 ^ z . (33) 

As far as only the case of small disturbance of the semiconfining regime of motion is 
of our concern, we can solve Eqs. (33) by perturbation theory writing the solutions 

as 

x = xi + x 2 , z = zi + z 2 , (34) 

where x\ and z\ are the solutions of the unperturbed equations (20) and x 2 and z 2 
are the solutions to the equations 

d 2 x 2 A o . . d z 2 <- 1 \s,~f2 /- in i-\ 

= -Su^x-i + x 1 z 2 ) +7^ , -p- = 2A Su x z x z 2 + 7^ . (35) 

From here, similarly to the way of Ref. [16], we get 



x 2 (t)= fG x {t-t') 
Jo 



7&(0 + ^Su; 2 x lZ2 (t') 



dt' , 



z 2 {t)= f ' G^t-t'^Ut')^ , (36) 
Jo 
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where 



G x (t) 



smh(pt) 



G z (t) 



smh((3 z t) 



p 



A 



1/2 



P Z = 2\\\P. 



Because of Eqs. (32), we have <C x 2 3> = <C z 2 ^> = 0. And for the mean square 
deviations, we obtain 



o 7 D x t 

< X 2 > = 



sinh(2/ft) 



- 1 



X 



smh(&t) [ cogh ^^^ _ cosh (^)] _ ^ S inh(/5 2 t)sinh(/5t) + cosh(&t)cosh(/%) - ll 
p z t p J 

(37) 

(38) 



2 7 D z t 
< z 2 > 



& 2 



sinh(2/3 z t) 



- 1 



These formulas generalize the result of Ref. [16]. However the main aim in con- 
sidering the role of random collisions is not just to derive formulas as above but 
rather to find conditions under which these collisions would not essentially disturb 
the semiconfined motion. Such a condition can be written as 



< 1 , D = sup{ A,, D y , D z } . 



(39) 



For estimates, we can take the collision rate as 7 ~ hpa s /m, where p is the density of 
atoms and a s , scattering length, and for the diffusion rate we may write D ~ k B T/h, 
where T is temperature. Then condition (39) yields 



(pai) 2 k B T 



2 \ 2 



mat 



(40) 



This shows that random atomic collisions will not disturb much the organized semi- 
confined motion if density and temperature are small enough to satisfy condition 
(40). 
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